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Non-Abelian Dual Superconductor Picture for Quark Confinement 
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We give a theoretical framework for defining and extracting non-Abelian magnetic monopoles 
in a gauge-invariant way in SU(N) Yang-Mills theory to study quark confinement. Then we give 
numerical evidences that the non-Abelian magnetic monopole defined in this way gives a dominant 
contribution to confinement of fundamental quarks in SU(3) Yang-Mills theory, which is in sharp 
contrast to the SU(2) case in which Abelian magnetic monopoles play the dominant role for quark 
confinement. 
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I. INTRODUCTION 

What is the mechanism for quark confinement? The 
dual superconductor picture proposed long ago [l| is be- 
lieved to be a promising mechanics for quark confine- 
ment. For this mechanism to work, however, magnetic 
monopoles and their condensation are indispensable to 
cause the dual Meissner effect leading to the linear po- 
tential between quark and antiquark, namely, area law of 
the Wilson loop average. 

The Abelian projection method proposed by 't Hooft 
0, Q can be used to introduce such magnetic monopoles 
into the pure Yang-Mills theory even without mat- 
ter fields, in sharp contrast to the 't Hooft-Polyakov 
magnetic monopole in the Georgi-Glashow gauge-Higgs 
model with adjoint matter fields. Indeed, numerical ev- 
idences supporting the dual superconductor picture re- 
sulting from such magnetic monopoles have been accu- 
mulated since 1990 in pure SU(2) Yang-Mills theory [J- 

However, the Abelian projection method explicitly 
breaks both the local gauge symmetry and the global color 
symmetry by partial gauge fixing from an original non- 
Abelian gauge group G = SU{N) to the maximal torus 
subgroup, H = U{1)^~^. Moreover, the Abelian dom- 
inance j4!| and magnetic monopole dominance [5| were 
observed only in a special class of gauges, e.g., the max- 
imally Abelian (MA) gauge [7] and Laplacian Abelian 
(LA) gauge, realizing the idea of Abelian projection. 

For G = SU{2), we have already succeeded to set- 
tle the issue of gauge (in) dependence by introducing a 
gauge-invariant magnetic monopole in a gauge indepen- 
dent way, based on another method: a non-Abelian 
Stokes theorem for the Wilson loop operator [^, 0] and 
a new reformulation of Yang-Mills theory rewritten in 
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terms of new field variables [Tol - IT2l | and [13l - IT5l |. elabo- 
rating the technique proposed by Cho [Iq and Duan and 
Ge [I4I independently, and later readdressed by Faddeev 
and Niemi [1^ and Shabanov |19| . 

For G — SU{N), N > 3, there are no inevitable rea- 
sons why degrees of freedom associated with the maximal 
torus subgroup should be most dominant for quark con- 
finement. In this case, the problem is not settled yet. 

In this paper, we give a theoretical framework for de- 
scribing non-Abelian dual superconductivity, i.e., super- 
conductivity caused by non-Abelian magnetic monopoles 
in D-dimensional SU{N) Yang-Mills theory, which 
should be compared with the conventional Abelian 
dual superconductivity, i.e., superconductivity caused by 
C/(l)^^^ magnetic monopoles in SU{N) Yang-Mills the- 
ory, hypothesized by Abelian projection. We demon- 
strate that an effective low- energy description for quarks 
in the fundamental representation can be given by a set 
of non-Abelian restricted field variables and that non- 
Abelian U{N — 1) magnetic monopoles in the sense of 
Goddard-Nuyts-Olive and Weinberg [2^ are the most 
dominant topological configurations for quark confine- 
ment as conjectured in 21, 22]. 

This paper is organized as follows. 

In section II, we rewrite the SU{N) Wilson loop op- 
erator in terms of a pair of gauge-invariant magnetic- 
monopole current k {{D ~ 3)-form) and the associated 
geometric object Ss defined from the Wilson surface E 
bounding the Wilson loop C, and another pair of an elec- 
tric current j (one-form independently of D) and the 
associated topological object Ny,, which follows from a 
non-Abelian Stokes theorem for the Wilson loop opera- 
tor [H]. 

In section III, we reformulate the SU{N) Yang-Mills 
theory in terms of new field variables obtained by change 
of variables from the original Yang-Mills gauge field 
■^1^ {x) J4J , so that it gives an optimal description for the 
non-Abelian magnetic monopole defined from the SU{N) 
Wilson loop operator in the fundamental representation 
of quarks. 

In section IV, we construct a lattice version of the refor- 
mulated Yang-Mills theory [23, l2g| to perform numerical 



simulations. The results of numerical simulations of the 
lattice SU{3) Yang-Mills theory give numerical evidences 
that the restricted field variables become dominant in the 
infrared for correlation functions and the string tension 
{infrared restricted non-Ahelian dominance) and that the 
U{2) magnetic monopole gives a most dominant contri- 
bution to the string tension obtained from SU{Z) Wil- 
son loop average [non-Ahelian magnetic monopole dom- 
inance) . 

The final section is devoted to conclusion and discus- 
sion. We will mention a possible direction of future works 
needed to confirm the non-Abelian dual superconductiv- 
ity. 



II. A NON-ABELIAN STOKES THEOREM FOR 
THE WILSON LOOP OPERATOR 

A version of a non-Abelian Stokes theorem (NAST) for 
the Wilson loop operator originally invented by Diakonov 
and Petrov [s!] for G = SU{2) was proved to hold 3 and 
was extended to G = SU{N) [2l|, [231 '^^ ^ unified way 
[23| as a path-integral representation by making use of 
a coherent state for the Lie group. For the Lie algebra 
su(A^)-valued Yang- Mills field £/f,{x) = s^i^{x)Ta with 
su{N) generators Ta {A ^ I,- ■ ■ ,N'^ - 1), the Wilson 
loop operator is defined by 



WcW] :=tr 



^ exp <^ igYM f 



dx'^Ji/f.ix) 



/tr(l). (1) 



The NAST enables one to rewrite a non-Abelian Wilson 
loop operator Wc \si\ in terms of an Abelian-like poten- 



tial Ay, as 



VFcM = / dnc{g)e^p 



iffYM 



A 



c 



(2) 



where gym is the Yang-Mills coupling constant, 
d/Jcig) '■= Yixec '^l^idx) '^ith s-n invariant measure d/i 
on G normalized as / d^(gx) = 1, g^; is an element of a 
gauge group G (more precisely, a representation Dfi{gx) 
of G), and the one-form A :— Ay{x)dx'-'' is defined by 

Ay{x) = tr{p[gls!/y{x)gx + igyugldygx]}, gx e G. (3) 

Here p is defined as p :— |A)(A| using a reference state 
(highest or lowest weight state of the representation) |A) 
by the use of a representation of the Wilson loop we con- 
sider. Note that tr(/c)) = (A|A) = 1 follows from the nor- 
malization of I A). Then it is rewritten into the surface- 
integral form using a usual Stokes theorem: 



W^cK] = / dp^{g)exp igyu / F 



where dp^{g) ;= Ilxes'^MCff^)' ^^"^ two-form F := dA 
\Fy^{x)dx^^ A dx'^ is defined by 



(4) 



Fy,{x) = ^2{N - l)/N[%,ix) + tgYlMpgl[dy,d,]gx}], 

(5) 



with the field strength '^^i, defined by 
%,^{x) ■■= dytT{n{x)£/^{x)} - d^tT{n{x)£/y{x)} 

+ j^ igyl^ti{n{x)[df,n{x),d^n{x)]}, (6) 

and a normalized traceless field n{x) called the color field 



n{x) := ^N/[2{N~l)]gx [p ~ l/tr(l)] gl (7) 

Finally, the Wilson loop operator in the fundamental rep- 
resentation of SU{N) is cast into the form [2J|: 



Wc 



dp^{g) exp {igyuik, ^s) + «.9ym(j, ^s)} , 



k:=5*f = *df, j:=Sf, f:^^2{N-l)/N^, 
Se := ^de^A-' = 5 * GsA-i, N^ := 6e^A-\ (8) 

where two conserved currents, "magnetic-monopole cur- 
rent" k and "electric current" j, are introduced, A := 
dS -\- 6d is the ZJ-dimensional Laplacian in the D- 
dimensional Euclidean space, and Q is an antisymmetric 
tensor of rank two called the vorticity tensor: Q'^^(x) := 
J^d'^S^'^{x{(T))S^{x — x{a)), which has the support on 
the surface S (with the surface element dS'^" {x{a))) 
whose boundary is the loop G. Incidentally, the last 
part iq yj.A .i^ p gj \dn. , d^] q.r I in F corresponds to the Dirac 
string 23, l2g| , which is not gauge invariant and does not 
contribute to the Wilson loop in the end. 

For SU{3) in the fundamental representation, the 
lowest-weight state (A| = (0,0, 1) leads to 

n{x) = gx{X8/2)gi G SU{3)/[SU{2) x f/(l)] ^ CP^, 

(9) 
with the Gell-Mann matrix As := diag.(l, 1, — 2)/V3, 
while for SU{2), (A| = (0, 1) yields 

n{x) = gx{'J3/2)gl e 5t/(2)/C/(l) c, S^ c, CP\ (10) 

with the Pauli matrix a^ := diag.(l, —1). The existence 
of magnetic monopole can be seen by a nontrivial Homo- 
topy class of the map n from the sphere S'^ to the target 
space of the color field n [2l|: For SU{3), 

n2{SUi3)/[SU{2) X t/(l)]) = ni{SU{2) x t/(l)) 
= 7ri(C/(l))=Z, (11) 



while for SU{2) 

7T2{SU{2)/U{1)) ^ 7Tl{U{l)) ^ Z. 



(12) 



For SU{3), the magnetic charge of the non-Abelian mag- 
netic monopole obeys the quantization condition |23j : 



Qm '■= / d^xk^ = 2'K\/igYlin, n G Z. 



(13) 



The NAST shows that the SU{3) Wilson loop operator in 
the fundamental representation detects the inherent U{2) 
magnetic monopole which is SU{3) gauge invariant, see 



(120)1 . The representation can be classified by its stability 
group H oi G [2l|, [23| • For the fundamental representa- 
tion of SU{3), the stability group is C/(2). Therefore, the 
non-Abelian U{2) ~ SU{2) x U{1) magnetic monopole 
follows from the field in the representation with the sta- 
bility group H = SU{2)i^2,3 X U{l)s, while the Abelian 
U{1) X U(l) magnetic monopole comes from that with 
H ~ U{1)3 X U{l)s- The adjoint representation belongs 
to the latter case. The former case occurs only when 
the weight vector of the representation is orthogonal to 
some of root vectors. The fundamental representation 
is indeed this case. For SU{2), such a difference does 
not exist and U{1) magnetic monopoles appear irrespec- 
tive of the representation, since H is always U{1) for any 
representation. For SU{3), our result is different from 
Abelian projection in which two independent U{1) mag- 
netic monopoles appear for any representation, since 

TT2iSU{3)/U{l) X C/(l)) = TTl{U{l) X [/(!)) = Z2. (14) 



III. REFORMULATING THE YANG-MILLS 
THEORY USING NEW VARIABLES 

Recently we have proposed a new reformulation [2J| 
of the SU{N) Yang-MiUs (YM) theory based on new 
variables by extending the Cho-Faddeev-Niemi (CFN) 
decomposition for TV > 3 [2^, [30|- Our reformula- 
tion allows options discriminated by the stability group 
H of the gauge group G — SU{N). When H agrees 
with the maximal torus group H = [/(1)^~^, it repro- 
duces a manifestly gauge-independent reformulation of 
the Abelian projection represented by the well-known 
maximal Abelian gauge. This case is called the maxi- 
mal option and agrees with the conventional CFNS de- 
composition for the SU{N) Yang-Mills theory for A^ > 3 
[29l ISOj . It was found ^2^] that there arc the other options 
with the stability group H other than the maximal torus 
group H = U{1)^~^. Such possibilities are overlooked 
so far. Especially, the case oi H = U{N — 1) is called the 
minimal option. The minimal option gives the optimal 
description of quark in the fundamental representation 
combined with the non- Abelian Stokes theorem for the 
Wilson loop operator given above. 

The reformulation enables one to understand quark 
confinement based on the dual superconductivity picture 
in a gauge independent way. This is because we can 
define gauge- invariant magnetic monopoles which are in- 
herent in the Wilson loop operator. 

In the realistic case of SU{3), there are two options: 
the minimal option [24l| with a single type of non-Abelian 
magnetic monopole characterized by the maximal stabil- 
ity subgroup ^ = U{2) = SU{2) x C/(l), and the max- 
imal one [2^, [so] with two types of Abelian magnetic 
monopoles characterized by the maximal torus subgroup 
H^U{1) X t/(l). 

We consider the decomposition of £/^{x) into two 



pieces yi_i{x) and ^^{x): 

<(a;) = r,,{x) + ^,,{x), 



(15) 



where the decomposed fields 7^ (a;) and ^(x) are deter- 
mined by solving the following defining equations (I) and 
(II), once a color field n{x) is given. 

(I) n{x) is a covariant constant in the background ^^{x): 

= D^[y]n{x) := d^n{x) - igYM[Mx),n{x)], (16) 

(II) ^^{x) does not have the i7- commutative part: 

^'{^)h ■■= (l - 2^^[", [n, •]]) ^'-(x) - 0. (17) 

From the viewpoint of quark confinement, this decom- 
position has a remarkable property: It is shown [3l[ that 
(II) guarantees (a), while (I) guarantees (b). 
(a) Yf^ alone reproduces the Wilson loop operator: 



Wc[£/]=Wc[ri (18) 



^cxp zgYM$ dx^r^ix)} /tr(l), (19) 



where 

Wc[y] :=tr 



(b) the field strength ^^^[T] := 8,,% - 9,,^^ - 
*5ym[^, "^1 in the color direction n agrees with 5f^i,: 



f^,{x)^tY{n{x).^^,Ar]{x)}. 



(20) 



This fact (a),(b) is also checked by using the explicit form 
of decomposed fields which are uniquely fixed by solving 
the defining equation: 



^^{x) ^ - tgyl ^^^^ ^^ [n{x), %[s^]n{x)] e £^e(G/^), 

^^{x) -.^^.{x) ~ '^^^~^\ n{x), [n{x),^^{x)]] € Cze{H), 
,^^{x) ■.=ig^l, ^^^'^\ n{x),d^n(x)] e Cie{G/H). 



N 



(21) 



In what follows, Cie{G) denotes the Lie-algebra of the 
Lie group G. ^ 



^ If the color field Ti(a;) is fixed to be a constant (this is considered 
to be a gauge fixing), then we have ^^(x) = 0. Then the decom- 
position ^fi{x) = ^(a;)-|-^p(a;) is reduced to the orthogonal de- 
composition: f^{x) = 'iffj.{x) = .b/^{x) ^ ~ ' [n, [n,.s^^(x)]] G 

Cie{H) and 3^^(x) = ^^^^^[n, [n,£i^^,{x)\\ G Cie{G/H). In the 
5(7(2) case and the maximal option of the SU{N) case, this pro- 
cedure under the reduction condition II25I I is equivalent to taking 
the Maximal Abelian gauge in the original Yang-Mills theory. 
But, the resulting theory in the minimal option of SU{N) case 
has never been considered to the best of author's knowledge. 
This theory will be studied elsewhere. 



By combining (a) and (b) with the NAST given in the 
previous section, therefore, the Wilson foop operator can 
be rewritten in terms of new variables: 



Wc[ 



dfis{g)exp igYMy/'2iN-l)/N 



S:dT,=C 



tr{n^[r]} . (22) 



In our reformulation, '^^^{x) and ^^{x) must be ex- 
pressed in terms of J2/^(a;). Therefore, we must give a 
procedure of determining n from ^, thereby, all the new 
variables "^p, <^ and n are obtained from ^ through 
the transformation law (|2T]): 



(23) 



Inserting 1 = ^ 'Dn°'Y[^5{n{x) - g^{\s/2)gl.) yields 
Wc ^Zyl, J dficig) J T^^^^T^n''S{n{x) ~ 9.{\»/2)gl) 

y g-S'YM[i;«']giSYM ic '^. (28) 

Thus, in the reformulated theory in which n^{x), '^^{x), 
Sy^{x) are independent field variables, Wc is written 



X p-SYM[n.'^,3'.'] igYMik,Sj:)+igYM{j,Ns) 



JgYMik,Sj:)+igYM{j,Ns) 



.//^(n^-^^jr,"). 



YM, 



(29) 



To solve this issue, we begin with counting degrees of free- 
dom: ^p e Cie{G) = su{N) means #[j/^] = D-dimG = 
D{N^ - 1), 'rf^ e Cie{H) = u{N - 1) means #[' 



where the Yang-Mills action is rewritten in terms of new 
variables using ([T5|) and ([2T|). Syuin,'^, ^] = S'ym[-2^] 
and the new partition function is introduced: Zym = 
fV'rf^V,T^^Vn^S{x)A'^^Je-^™^'^''^'-^l It is shown 



means #[«''] = dim(G/7J) = 2(A^ - 1). Thus, the new 
variables (nP,^!;, JTj") have the 2{N — 1) extra degrees 
of freedom, to be eliminated to obtain the new theory 
equipollent to the original one. For this purpose, we im- 
pose 2{N — 1) constraints X = 0, which we call the re- 
duction condition. For example, minimize the functional 



R\. 



j d^x\{D,[s^]nf 



(24) 



D ■ dimff = D{N - 1)2 and JT^ € Cie{G / H) means \m that the integration measure Vs/^ is finally trans- 
#[S'^f'] = D.dim(G/i7) = 2D{N-l) and n e Cie{G/H) formed to Vt^^V3C^VnH{x)Af^J, where (i) the Jaco- 

bian J is very simple, J = 1, [23| irrespective of the choice 
of reduction condition, (ii) ;\;[^,n] = is rewritten in 
terms of new variables: x := x["/^, ^] := D^'[r]3^^, 
and (iii) the associated Faddeev-Popov determinant App 
is calculable using the BRST method, e.g. ^11)]. 

In the previous section, the Wilson loop operator has 
been exactly rewritten in terms of the gauge-invariant 
magnetic current k (and the electric current j). This 
shows that the Wilson loop operator can be regarded as a 
probe of magnetic monopoles. In this section, moreover, 
we have succeeded to connect the Wilson loop average 
with magnetic monopoles which is supposed to be a basic 
ingredient to cause dual superconductivity as a promis- 
ing mechanism of quark confinement. In fact, Eq. (j29p 
tells us what quantity we should examine to see the mag- 
netic monopole contribution to the Wilson loop average. 
This equation is important to give a connection between 
our formulation and magnetic monopole inherent in the 
Wilson loop to see quark confinement. In fact, we give 
numerical calculations of the potential Vm (R) in the next 
section based on the lattice version of Eq. 



with respect to the enlarged gauge transformation: 
Df^[£/]u), and Sn = gi[0,n] — gi[9±,n] where 



Ss^„ 



U3 e Lie{G) and 6»^ G Cie{G/H). Then, we find 
5R[s^,n\ = gjd^x{ei_ - u±) ■ i[n,Dt'[s^]Df,[£/]n], 
where tjj± denotes the component of u; in the direction 
£ie{G/H). The minimization dR[£/,n] = imposes no 
condition for u)± = 6± (diagonal part of G x G/H), while 



x[-^,n] :- [n,D''[,^]D^[^]n] - 0, 



(25) 



is imposed for u}± ^ 0± (off-diagonal part of G x G/H). 
The number of constraint is #[x] = dim(G x G/H) — 
dim(G) = dim(G/i^) = 2{N - 1) = #[h^] as desired. 
As a bonus, the color field n{x) is determined by solving 
(^5)1 for given ^(x). This completes the procedure. 
The Wilson loop average Wc is defined by 



Wc = (VFcH)ym := ^YM 



A -SymIs^] 



WcW], 
(26) 



with the partition function Zym = /X>i2^/e~-^YM[.s/] ^y 
omitting the gauge fixing to simplify the expression. The 
pre-NAST ^ tells us that 

Wc = ^YM / rfMc(5)2?^/e-^™l-^le*9™ fc ^. (27) 



IV. LATTICE SU(N) YANG-MILLS THEORY 

A. Reformulating the lattice SU{N) Yang-Mills 
theory 

We have reformulated the Yang-Mills theory with a 
gauge group G = SU{N) also on a lattice using new 
variables in the same spirit as in the continuum, see |25l . 
l26j for the details. But we just summarize it below. 

First, gauge field configurations {Ux,fi} on a four- 
dimensional Euclidean lattice are generated by using the 
standard method: the Wilson action and pseudo heat- 
bath algorithm. The gauge variable Ux,fi on the link has 



the gauge transformation: 



'~'x,^i ~^ ^''x'^x,fj,''''x+fi ^ ^X.fJ.^ 



n^ e G. 



(30) 



Second, according to the continuum formulation [2J], 
we introduce just a single color field n^ even for G = 
SU{N) {N > 2) in the minimal option. The color field 
Ux on a lattice is regarded as a site variable defined on a 
site X and takes the value as 



rix = n^TA e G/ii ^ SU{N)/U{N - 1) ~ CP 
with a unit length 

n^n^ = 1. 



iV-l 



(31) 



(32) 



For a given set of gauge field configurations {Ux,fi}, a set 
of color fields {Ux} is determined by imposing a lattice 
version of the reduction condition. A reduction condition 
in the minimal option on a lattice is given by minimizing 
the reduction functional Prod for a given set of gauge 
field configurations {Ux,^} with respect to the color field 
{nx}- 

F,ed[n,C/]:-e^^tr{(i?^[C/]n,)(i?^[C/]n,)t}/tr(l) 

X,fJ, 

=£^-2 ^[^ _ 2tr{nxUx,^nx+^Ul^)], (33) 



x,// 



where Dl^[U] is the lattice covariant derivative in 
the adjoint representation defined by D^^[U]nx := 



~\Ux„ 



,n 



x+pi 



nxUx^fj.) with a lattice spacing e. Thus, 



a set of color fields n{x) we need is obtained as a set of 
unit vector fields n(x) which realizes the minimum of the 
reduction functional: 



Prcd [n,U]^ min F^d [n,U], 



(34) 



It is observed that solving the reduction problem is equiv- 
alent to finding the ground state of the spin-glass model, 
since the reduction functional F,-cd is minimized with 
respect to the color field {rix} under the random link 
interaction Jx^[U] for given gauge field configurations 
{Ux,f.}: 



F,ed[n,t/] ^e^-'Y.(^-J^,^[U]nS^, 



J^^^iU] ■.^2tr{TAUx,f.TBUl^). 



(35) 



This observation has been actually used to find the min- 
imum in SU(2) case [32] ■ After applying the reduction 
condition [25| , the color field transforms under the gauge 
transformation in the adjoint way: 



Ux -> ftxnxftx ^ = "L: ^x e G. 



(36) 



The reduction functional F^^d is invariant under the 
gauge transformation. Therefore, imposing the reduction 
condition does not break the original gauge invariance. 



We can impose any gauge fixing afterwards, if necessary. 
The details for the algorithm of the reduction jarocedure 
on a lattice in the SU(3) case will be given in [33] . 

Third, new variables on a lattice are introduce d by 
using the lattice version of change of variables [2^ [2a] : 
Once a set of color fields rix is given, the G = SU{N)- 
valued gauge variable Ux,p, G G is decomposed into the 
product of two G- valued variables Xx^n and Vx^/^ defined 
on the same lattice: 



Ux.fi — Xx,fj.Vx,fi G G, Xx^fj,, Vx^f^ G G, 



(37) 



where the lattice variables Vx,^ and Xx,f_i are supposed 
to be related to the Lie-algebra ^(a;) and ^fj,{x) as 

Vx,f, = exp{-iegYM'>^/x(a;)}, Xx,f_c = exp{-iegYM^t^{x)}, 

(38) 
just as 



Ux.f, = exp{-ie5YM^/i(a:)}. 



(39) 



We require that T4^^ is a new link variable which trans- 
forms like a usual gauge variable Ux.^ on the same link: 

Vx,^^nxVx,^n-l^^v^^^, nxeG. (40) 

For this gauge transformation to be consistent with the 
decomposition ([37]), consequently, Xx^^ must behave like 
an adjoint matter field defined at the site x under the 
gauge transformation: 

Xx,^ ^ nxXx^^^n-' = x'^^^, nx e g. (4i) 

These properties of the decomposed variables under the 
gauge transformation are expected from the continuum 
version. The decomposed variables Xx.f^ and Vx^^ are 
determined by solving defining equations. A lattice ver- 
sion of the first defining equation proposed in |25l . 1261 ] 
is: (I) The color field Ux is covariantly constant in the 
(matrix) background Vx^fj,'- 



= £D'^^\V]nx ■■= Vx^^^rix+ti - rixVx., 



(42) 



where D)^' \V] is the lattice covariant derivative in the ad- 
joint representation. The solution of this defining equa- 
tion can be obtained exactly for any N 26] (without 
using the ansatz employed in [25| to find the solution for 
A^ = 2, 3) to give A^,,^ and T4,m = XljJx.^,: 



Jv-r 



--Ll(AetCLx,^)f"'9x\ 



Vx,t, =9xLx,t,Ux^^{Aei{Lx^^,)) ^/^, 



(43a) 
(43b) 



where 



-^x.,fi — [y -'-^x,i.i-^x,fj,) J-^x 



X,fl X,fl 



T 7-1" 

^X,^-'-'X.l.i 



(43c) 



with 



monopole current k defined in 



' N " 



1 



+ {N~2) \/ ^^^V^ K + U:,,^ri,+^U-l) 



N 



+ 4 (iV - 1) n^U^^^n.^+f,U^l. 



(43d) 



Here a common factor g^; in the above expressions for 
Xx^^ and Vx,^ is the part undetermined from the first 
defining equation alone. In fact, g^ is an element of the 
extra symmetry associated with the decomposition '2(^: 
Z{N) X H, H = U{N - 1) C SU{Ny. 



gx ^e-^'^*"-/^ exp < -zfl: 



(g, = 0,---,iV-l), 






(44) 



where ax, ax S K and {ux } is a set of Hermitian trace- 
less generators of SU{N — 1) commutable with Ux- 

In order to fix it, we must impose further conditions. 
Hence we impose the second defining equation, e.g., 
(II) gx is equated with an element gj!: 



9x =9: 



,0 



(45) 



The simplest one is to take g'^ — 1. Thus the decomposed 
variables Xx,f_i and Vx^fi are completely determined. It 
can be checked that the lattice formulation given in this 
section reproduces the continuum formulation given in 
the previous section, in the naive continuum limit e — > 0. 



B. Numerical simulations for SU{?i) Yang-Mills 
theory 

In the new formulation, we can define another non- 
Abelian Wilson loop operator Wc [^] by replacing s^ by 
y in the original definition of the Wilson loop opera- 
tor Wc[-s/]- For the lattice version of the Wilson loop 
operator Wc[£^], 



Wc[U]:=tr 



II Ux,x+fj. 



Ml), (46) 



the lattice version of Wc\y\ is easily constructed: 



Wc[l^]:=tr 



W Vx,x+fj. 

<x,x+fj.>^C 



Ml). (47) 



This is invariant under the gauge transformation ()40p . 
Moreover, we define the lattice version K of the magnetic 



2 



'■a:, /J. ■ — ^u ^^x^^y — cj ^/ii/ap 



(?Qx,ai3 :=arg tri f -1 

x\4,„v;+„,^i;V^,„<4]. (48) 



Vl"^ 



It is easy to observe that Qx,^ji.v is invariant under the 
gauge transformation ([55)) and (PH)) . and hence Kx,^ is 
also gauge- invariant. Then we can define the magnetic- 
monopole part of the Wilson loop operator by 

Wc[K]■.^e^^p(^Y,Kx.,^^^lA, 

\ x,fi ) 

S^,M -E^l'('^ - s')ie^a/j^a„5/^(s' + ii), (49) 



where Si As' + /i) is a plaquette variable satisfying 
dg,S^p{x) = Jp{x) with the external source J^,^ intro- 
duced to calculate the static potential, d' denotes the 
backward lattice derivative d fx — fx ~ fx-fj., S^ a de- 
notes a surface bounded by the closed loop C on which 
the electric source Jx.fj. has its support, and A^ {x — x') 
is the inverse Lattice Laplacian, see e.g., |35| . 

Numerical simulations are performed for SU(3) Yang- 
Mills theory on the 24^ lattice according to the lattice 
reformulation explained above for A^ = 3. More details 
of numerical simulations will be given in a subsequent 
paper |33j . 

The static quark-antiquark potential Vf (R) is defined 
by taking the limit T — > oo from the Wilson loop average 
(Wc[i7]) for a rectangular loop C — Rx T: 



Vf{R) 



hm -\n{Wc[U]). 

1 — ^OO 1 



(50) 



In practice |40|, we fit numerical data of (VFc[C/]) by the 
two- variable function W{R,T) according to 

{Wc[U])^eM-W{R,T)), (51a) 
W{R, T) -.^TViR) + (aiR + bi+ ci/R) 

+ {a2R + b2+C2/R)T-\ (51b) 

V{R):=aR + b + c/R, (51c) 

and determine all coefficients in W{R, T). Then we iden- 
tify Vf{R) with V{R) to be obtained by extrapolating 
W (R, T)/T to T^oo: 

V{R)=\un ^^^^Vf{R). (52) 

Here the coefficient a of the linear part of the potential 
(|51c[) is the string tension which equals to the slope of 
the curve for large R. 
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FIG. 2: Color field correlators (n'^(a;)n^(0)) {A,B = 
I,--- ,8) as functions of the distance r := |a;| measured at 
/? = 6.2 on 24* lattice, using 500 configurations under the 
Landau gauge. (Left panel) A = B, (Right panel) A ^ B. 



FIG. 1: 5'L'^(3) quark-antiquark potentials as functions of the 
quark-antiquark distance R: (from above to below) (i) full 
potential Vf{R) (red curve), (ii) restricted part Vr{R) (green 
curve) and (iii) magnetic-monopole part Vm,{R) (blue curve), 
measured at /3 = 6.0 on 24'' using 500 configurations where e 
is the lattice spacing. 



In Fig. [TJ we compare the three quark-antiquark po- 
tentials (i), (ii) and (iii). For each potential, we plot a set 
of point data for a specified value of T (e.g., T = 6, 10): 



— ln{Wc[-]) versus R, 



(53) 



and the curve 

ViR)=(TR + b + c/R, (54) 

extrapolated to T — ?► oo according to ([?T|) and ((5^ : 

(i) the full potential Vf{R) calculated from the stan- 
dard SU{3) Wilson loop average {Wc[U]), 

(ii) the restricted potential Vr{R) calculated from the 
decomposed variable f through the restricted Wil- 
son loop average (W^c[^]) 



Vr{R) = - \im hn{Wc[V]}, 

T— i-oo 1 



(55) 



(iii) magnetic-monopole contribution V,n (R) calculated 
from the lattice counterpart (|49|) of the continuum 
quantity (e'C"^"^)) according to ^^: 

V,n{R)=- hm ^\n{Wc[K]), (56) 



Three potentials are gauge invariant quantities by con- 
struction. 

The results of our numerical simulations exhibit in- 
frared restricted variable f' dominance in the string ten- 
sion, e.g.. 



av _ 0.0380 
ov ~ 0.0413 



0.92, 



(57) 



and non-Abelian magnetic monopole dominance in the 
string tension, e.g., 






0.0352 
0.0413 



0.85. 



(58) 



However, we know that ajj has the largest errors among 
three string tensions. Incidentally, if we use the other 
data for ey/aipF at /3 = 6.0 given in Table 4 of [34| where 
e^au* = (eV^V^)^ = 0.2154^ - 0.2209^ = 0.0464 - 
0.0488, the ratios of two string tensions av,o'M to the 
total string tension au are modified 



—^ ^ 0.78 - 


- 0.82, 


(59) 


^M ^Q^72^ 


- 0.76, 


(60) 



au* 

Anyway, we have obtained the infrared restricted vari- 
able 'f^ dominance in the string tension (78-82%) and 
the non-Abelian magnetic monopole dominance in the 
string tension (72-76%). Both dominance are obtained 
in the gauge independent way. ^ 

To obtain correlation functions of field variables, we 
need to fix the gauge and we have adopted the Landau 
gauge for the original Yang-Mills field ^ so that the 
global color symmetry is not broken. This property is 
desirable to study color confinement, but it is lost in the 
MA gauge. 

FiglH shows two-point correlation functions of color 
field {n^{x)n^{Q)) versus the distance r := |a;|. All plots 
of correlators for A = S = 1, 2, • • • , 8 overlap on top of 
each other, and hence they can be fitted by a common 
non- vanishing function D{r) (left panel), while all corre- 
lators for A ^ B are nearly equal to zero (right panel). 
Therefore, the correlators {n"^ {x)n^ {Q)) are of the form: 

{n^{x)n^{Q))^5^^D{r) (A, B = 1, 2, • • • , 8). (61) 

We have also checked that one-point functions vanish. 



Hx)) = ±0.002 ~ (A==l,2,- 



(62) 



^ The method of fitting the data given in this paper is the same 
as that in |40||, but is different from that used in [3al. 
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FIG. 3: Field correlators as functions of the distance r :— \x\ 
(from above to below) {f^f {x)^^ (0)) , {.<i/f{x)a/^f{0)}, and 
{^/(a;)^/(0)>. 



These results indicate the global SU{3) color symme- 
try preservation^ i.e., no specific direction in color space. 
This is expected, since the Yang-Mills theory should re- 
spect the global gauge symmetry, i.e., color symmetry, 
even after imposing the Landau gauge. 

Fig. 13] shows correlators of new fields y,^, and origi- 
nal fields £if . This result indicates the infrared dominance 
of restricted correlation functions ('^'^(a;)'^'^(0)) in the 
sense that the correlator of the variable Y behaves just 
like the correlator (j2/^(x)j2k^(0)) of the original variable 
£/ and dominates in the long distance, while the correla- 
tor ( jr'/(x)J"/(0)) oiSU{3)/U{2) variable jT decreases 
quickly in the distance r. For <^, at least, we can intro- 
duce a gauge-invariant mass term 



iMi^/j-^, 



(63) 



since ^ transforms like an adjoint matter field under 
the gauge transformation. In view of this fact, we fit the 
data of the contracted correlator (^'^(a;),^'^(0)) using 
the "massive" propagator for large r :— \x\: 

(jr/(x) jr-4(o)) = f£^e^'^ ^ 



(27r)4 ^ B + Ml 



.-Mxr 



const. 



j,3/2 " 

Then the naively estimated "mass" Mx of ^ is 



Mx = 2.409^/0^ = l.lGeV. 



(64) 



(65) 



This value should be compared with the result in MA 
gauge [36|. For more preliminary results of numerical 
simulations, see [39| for various properties of magnetic 
monopoles of SU(2), [33] for the maximal option of SU(3) 
and [3a, |40| for the minimal one of SU(3). 



V. CONCLUSION AND DISCUSSION 

In this paper, we have given a new reformulation called 
the minimal option of the SU{N) Yang-Mills theory in 
the continuum and on the lattice by elaborating previous 
works to provide one with an efficient framework to study 
quark confinement originating from magnetic monopoles 
defined in the gauge-independent manner. 

In fact, we have given first numerical evidences that 
non-Abelian magnetic monopoles defined in a gauge- 
invariant way in this paper are dominant for confine- 
ment of fundamental quarks in SU(3) Yang-Mills the- 
ory. By using the gauge invariant magnetic current fc, we 
have extracted just the U{1) part of the stability group 
U{N - 1) ~ SU{N - 1) X [/(I) for the non-Abelian mag- 
netic monopole associated with quarks in the fundamen- 
tal representation, which is consistent with the consid- 
eration of the Homotopy group (fTTj) . This [/(I) part is 
enough to extract the dominant part of the Wilson loop 
average. 

This fact suggests that the non-Abelian dual supercon- 
ductivity caused by condensation of non-Abelian mag- 
netic monopoles could be a mechanism for quark con- 
finement in SU(3) Yang-Mills theory. However, in or- 
der to establish the non-Abelian dual superconductiv- 
ity in SU(3) Yang-Mills theory in this sense, we must 
confirm that the spontaneous breaking of the dual U{1) 
symmetry is associated with condensation of the non- 
Abelian magnetic monopoles obtained in this paper. For 
this purpose, we need to specify how the relevant low- 
energy effective theory of the SU{3) Yang-Mills theory 
looks like, although it is believed to be a dual Ginzburg- 
Landau model based on the dual superconductor picture 
for quark confinement. 

In order to see directly the non-Abelian nature of mag- 
netic monopoles defined in this paper, it will be neces- 
sary to study how the interaction among the non-Abelian 
magnetic monopoles is described in the short distance 
where the internal non-Abelian SU{2) degrees of freedom 
in f/(2) ~ SU{2) X C/(l) other than [/(I) are expected to 
become relevant. Moreover, it will be interesting to study 
the contribution of non-Abelian magnetic monopoles to 
gluon confinement. In the framework of the reformula- 
tion using new variables, a first step in this direction was 
taken recently in [4i| for SU(2) gauge group. Extend- 
ing this work to the SU{3) gauge group will be a main 
subject of subsequent works. 
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